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ABSTRACT 

We present both phenomenological and more physical photoionization models of the Chandra HETG 
spectra of the Seyfert-1 AGN NGC 4051. We detect 40 absorption and emission lines, encompassing 
highly ionized charge states from O, Ne, Mg, Si, S and the Fe L-shell and K-shell. Two indepen- 
dent photoionization packages, XSTAR and Cloudy, were both used to self-consistently model the 
continuum and line spectra. These fits detected three absorbing regions in this system with densities 
ranging from 10 10 to 10 11 cm~ 3 . In particular, our XSTAR models require three components that 
have ionization parameters of log£ = 4.5, 3.3, & 1.0, and are located within the BLR at 70, 300, 
and 13,000 R g , respectively, assuming a constant wind density. Larger radii are inferred for density 
profiles which decline with radius. The Cloudy models give a similar set of parameters with ionization 
parameters of log£ = 5.0, 3.6, & 2.2 located at 40, 200, and 3,300 R g . We demonstrate that these 
regions are out-flowing from the system, and carry a small fraction of material out of the system rel- 
ative to the implied mass accretion rate. The data suggest that magnetic fields may be an important 
driving mechanism. 



1. INTRODUCTION 

The advent of CCD spectro scopy with ASC A showed 
that the "warm absorbers" (jHalperrJ 1984) that are 
present in the X-ray spectra of some Seyfert galaxies can 
be described with a combination of He-like O V II and 
H-like O VIII absorption edges (jRevnoldsl fl997t) . The 
spectral resolution of CCDs is insufficient to clearly re- 
veal individual narrow absorption lines that should be 
associated with such edges. The high resolution disper- 
sive spectrometers aboard Chandra and XMM-Newton 
have been able to clearly reveal such narrow absorption 
lines in a number of Seyfert AGN. Indeed, it is now clear 
that X-ray warm absorbers span a range of ionizations, 
generally have blue-shifts on the order of 10 2 ~ 3 km s _1 
with respect to their host galaxy, and so must be outflows 
(|Kaastra et al.ll2000HKaspi et al.ll2000l I2002T) . 

With deep observations and careful, self-consistent 
photoionization studies of high resolution spectra, the 
nature and constraints on the origin of these outflows can 
- in principle - be determined. In some warm absorbers, 
there is evidence that the opacity must derive partl y from 
elements trapped on dust grains (ILee et al.ll200l . sug- 
gesting that the absor ption may partia lly originate in 
a parsec-scale "torus" ( A ntonuccl 119931 ) . Similarly dis- 
tant absorption is a via ble description of the absorbers in 
many spectra (see, e.g. [Blustin et al1l2005|) . These out- 
flows are plausibly powered by radiation pressure from 
the central engine. In this picture, the outflows may be 
termed "disk winds", but they are unlikely to provide 
deep insights about the nature of disk accretion for they 
lack the ability to transfer angular momentum out of the 
system. 

The best spectra of Seyfert warm absorbers may sug- 
gest a different origin. The 900 ksec Chandra/ 'HET l G 
spectrum of NGC 3783 is richly detailed, and photoion- 
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ization analysis suggests that it may originate within 
several hundred gravitational radii (r g = GM/c 2 ) of 
the black hole (jKaspi et al.l [2002}. Detailed spec- 
tral modeling of the X-ray absorption in NGC 4151 
finds high column densities and evidence of large bulk 
motion that may be require large non-radial motions 
(|Kraemer et al.l [2005f ). This is qualitatively consistent 
wit h magnetocentrifuga l disk winds of the kind described 
bv iBlandford fe Pavnel (|1982[ ). Owing to the fact that 
these winds can remove angular momentum from ac- 
creting gas (thus allowing it to accrete), some warm ab- 
sorbers (or the most ionized, innermost regions of warm 
absorbers) may hold clues to the fundamental physics 
of disk accretion (see, e.g. iProga fe Begelmanl |2003; 
iFukumura etail [2010) for newer theoretical models). 
Even stronger evidence of magnetic processes driving 
disk winds and a possible anti-correlation with radio jets 
has been fou nd in the observations of some stellar-mass 
black holes (iMiller et all l2006al 120081: iKallman et al.l 
l2lMlLTiketic et al.ll2010t Blum et al.ll2010D . 

Understanding the ultimate influence of these warm 
outflows on host envir o nmen ts is also interesting and 
timely. iBlustin et al.l (|2005[ ) considered a sample of 
Seyfert warm absorbers, and derived typical values for 
these outflows. If the winds operate during the entirety of 
the active Seyfert/ AGN phase, then warm absorbers may 
expel up to 10 8 M Q of material into their host galaxies. 
In large elliptical galaxies, this may be unimportant com- 
pared to the stellar mass. However, this may represent a 
substantial fraction of the gas in the central spheroid of 
spiral galaxies. 

Estimates of the radii at which warm absorbers exist, 
and their total mass outflow rates, depend on detailed 
spectroscopy. Factors such as the gas density, covering 
factor, and filling factor must be estimated using line 
ratios and a number of physical assumptions. Indeed, 
the geometry of warm absorber outflows is poorly un- 
derstood; at present, only lines in N GC 3783 show eve n 
modest evidence of P-Cygni profiles ( Kaspi ct al. 2002) , 
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likely indicating that the wind is quasi-spherical. It is 
notable, then, that the best spectra provide the best evi- 
dence that warm absorbers may be disk winds that could 
be related to fundamental disk physics. 

Clearly, deep observations and sensitive spectra are 
vital to understanding accretion onto, and mass out- 
flows from, the central engine in Seyfert AGN. To this 
end, we have analyzed the summed, deep-exposure Chan- 
dra/HETG spectrum of the Seyfert- 1 galaxy NGC 4051. 
This Seyfert is close (z — 0.00234), and the mass of 
the black hole has been constrained using reverbera- 
tion mapping t e chniq ues (M B h = 1-9 ± 0.8 x 10 6 M Q , 
Peterson et all (l200l) : M BH = 1.7 ± 0.5 x 10 6 M@, 
Dennev et al.l (|2009t )). The innermost accretion flow in 



NGC 405 1 has been studied using relativisti c disk reflec- 
tion (e.g.. lOgle et al.ir2004J : iPonti et al.ll2006l or a review, 
see Miller 2007). 
The X-ray warm ab sorber in t his AGN has pre- 
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200 41 iNucita et al l 
Lobban et al.l 120111 : 
We confirm the multiple ab- 
sorbing regions with a range of ionization parameters, 
1 < log£ < 5 described in these papers. However, we 
do not find the various high velocity compone n ts, i.e . 
> 2000 km s " 1 , cla imed b y ISteenbrugge et al.l (2009), 
lLobban et all ()2011h . and iPounds fe Vaughanl (|2011h . 
Furthermore, we present evidence for material arising 
closer to the central engine than previously claimed, i.e. 
< 13000i? 9 . 

Herein, we present phenomenological and self- 
consistent, physical, photoionization models for the X- 
ray warm absorber observed in the time-averaged HETG 
spectrum of NGC 4051. We find that the wind likely 
originates close to the black hole, within a few hundred 
gravitational radii. We also find modest evidence of P- 
Cygni profiles in the H-like charge states, suggesting that 
the hottest, innermost part of the wind is quasi-spherical. 
A cosmology of H = 70 km s _1 Mpc -1 , qO =0 and A 
= 0.73 was used throughout this analysis. 

2. METHODS 
2.1. Data Reduction 

We have analyzed twelve Chandra High Energy Grat- 
ing (HEG) and Medium Energy Grating (MEG) spec- 
tra of the Seyfert- 1 NGC 4051, between the dates of 
09 November 2008 (MJD 54779) and 30 November 2008 
(MJD 54800). The event files, redistribution matrices 
and auxiliary matrices were created via the CIA0 com- 
mands (version 4.2) tgextract, mkgrmf and dmarfadd 
respectively. After acquiring the data and response files, 
we co-added the twelve HEG and MEG observations with 
the HEAS0FT, version 6.9, procedures mathpha, addrmf , 
and addarf . Although NGC 4051 is know to be highly 
variable, especially in the X-ray band, these data did not 
vary in flux by more than a factor of two and in spectral 
index varies from 1.7 to 2.0. Finally, by utilizing both the 
positive and negative first order spectra, we were able to 
get over 0.17 ct/s and avoid pile-up typically associated 
with the zeroth order spectra. A total exposure time of 
308 ksec was achieved. 

We utilized the Interactive Spectral Interpretation Sys- 
tem (ISIS), version 1.6.1, as well as XSPEC, version 12.6.0, 
to analyze the time-average spectra employing both phe- 



nomenological and photoionization modeling techniques. 
The data were binned at a minimum of 10 counts per 
bin. 

2.2. Phenomenological Model 

Initially, we fit the HEG and MEG broadband spectra 
using Xspec between the ranges of 1.5-20A and 1.5-23A 
respectively. The ranges were chosen because the respec- 
tive sections had the highest sensitivity for each spectro- 
graph, ensuring a maximum signal to noise ratio in each 
spectrum. The preliminary model consisted of a power- 
law and a disk blackbody component. The disk black- 
body is not likely to be physical; however, it adequately 
describes the soft-exce ss seen in NGC 4051 and many 
Seyfert galaxies (e.g. , [Reynolds 1997; Magd ziarz et al.l 
119981: iCrummv et all 120061) . We also include an effec- 
tive H column density in all models with the component 
phabs, which is frozen at 1.15 xlO 20 cm" 2 to model the 
line-of-sight absorption (|Kalberla et al.l 12005). Finally, 
we included a multiplicative constant, so as to account 
for any discrepancy between the normalization of the two 
spectrographs. 

The best-fit parameters of the model with a \ 2 l v — 
9067/7531 are listed in Table HJ The ratio of the data to 
the continuum model is plotted in Figure Q] The power- 
law index of T = 1.94±0.01 is typical of Seyfert galaxies. 
The disk blackbody temper ature of 0.118± 0.005 keV is 
also typical of Scyfcrts fe.g.. lRevnoldslll997l ). but again, 
is not likely to be a physical model. 

In general, this continuum-only model is a poor fit to 
the data, and it is clearly due to the fact that emission 
and absorption features are present. We note them in 
the residuals in Figure [TJ We first employ a phenomeno- 
logical model as a means of characterizing both of these 
features. The spectra were fit with ISIS utilizing 2 or 
3 A sections with a local power-law continuum and a 
Galactic absorption of 1.15 xlO 20 atoms cm -2 . By al- 
lowing the continuum to vary from segment to segment, 
we were able to precisely detect any underlying features 
amongst the residuals. We used a Gaussian line profile 
to describe the residual absorption and emission features. 
A total of 40 lines were detected and are listed in Tables 
and [3] Of these lines, all but one are detected at least 
at the 3er confidence level, and the remaining line is de- 
tected at a 2cr confidence level, but it is expected as part 
of a pair. The significance of each detection, which is 
determined using Ftest statistics, is also listed in these 
tables. The column densities listed in Tables [3] are also 
derived for each ion with known oscillator strength using 
the relation W\ = (7re 2 /m e c 2 )Nj\ 2 fy, where W\ is the 
equivalent width, Nj is the column density, A is the rest 
wavelength, and is the oscillator strength. 

We detect a number of different, highly ionized charge 
states including species of O, Ne, Mg, Si, S and Fe L-shell 
transitions, as well as narrow fluorescent Fe Ka and f3 
lines and a broad Fe Ka line. The median velocity of all 
of these features is -580 km s _1 , while the mean velocity 
is -530 km s _1 with a standard deviation of 320 km s _1 . 
The emission features have a slightly lower mean velocity 
of -310 km s _1 with a standard deviation of 450 km s , 
consistent with emission in the rest frame of NGC 4051. 
The absorption features have a faster mean outflowing 
velocity of -630 km s" 1 with a standard deviation of 200 
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km s~ 4 . The median and mean measurements do not 
include the unidentified lines or the fluorescence emission 
lines. 

These velocity shifts suggest that the absorption fea- 
tures are produced in a region along our line-of-sight that 
is moving toward us. Furthermore, the fact that the ab- 
sorption is blue-shifted with respect to the emission fea- 
tures may suggest the evidence of P-Cygni lines. These 
features are most clearly seen in Lyman series of both 
O VIII and Ne X, in particular the lines with rest wave- 
lengths of 18.96 A, 16.01 A and 12.1 A (see Figures M 
and [3]). It should be noted that any emission features 
in the He-like forbidden or intercombination lines can- 
not be P-Cygni features, due to the fact that they are 
optically thin. However, the detection of an absorption 
feature blue-ward of the emission feature in the H-like 
lines may give insight to the geometry of this system in 
the inner highly ionized region. If these two features are 
connected and not coincidental, then the lines suggest 
a spherical shell of material that is out-flowing from the 
system, and lends itself to determining roughly the cover- 
ing fraction of this system. However, we cannot rule out 
that the H-like lines may also be consistent with emis- 
sion in the rest frame of NGC 4051 and an absorbing 
region further from the central engine. To date, the only 
evide nce for P-Cygni pr ofiles in Seyfert winds is in NGC 
3783 (jKaspi et alJl2002l) . 

2.3. Density Diagnostics 

There are four particularly relevant sets of line detec- 
tions found in this data. The first two are the He-like 
triplet emission lines of O VII and Ne IX, and remaining 
two are the Fe L-shell doublets, Fe XXI and Fe XXII. The 
line ratios in the triplets allow for temperature and den- 
sity diagnostics. However, this assumes that the emission 
lines originate in the same region as the absorption fea- 
tures and line ratios ar e not also affected by the ultravio- 
let radiation field (e.g.. lMauche et al. 2001$. In addition, 
in the spectrum of NGC 4051, the emission component of 
the resonance line (r) is badly blended with the absorp- 
tion component, as indicated by an absorption equivalent 
width of O VII 18.65A which is larger than the equiva- 
lent width of the 21. 6A line (Figured]). Therefore, the 
well-known G ratio is unobtainable, giving only the ratio 
of the forbidden (f) to intercombination (i) lines. The 
ratios 

R(n e ) = ! - (1) 
i 

are 3.4 ± 1.6 and 1.1 ± 0.4 for the O VII and Ne IX 
R ratios. These R ratios are both consistent with the 
R ratios usin g the O VII and N VI triplets i n work by 
iNucita et all (|201Q[ ). However, INucita et all (|2010l ) do 
find a lower limit to the ratio of R > 5.6 for the Ne IX 
triplets which is much higher then our results. 

The R ratios we present would imply densities of 
10 10 cm ~ 3 and 10 12 cm ~ 3 for the ratios from O VII 
an d Ne IX, respectively, ac cording to the calculations 
of iPorquet fe Dubaul ([2000) if the UV flux is not im- 
portant. However, photons near 1620 A and 1248 A 
can excite transitions from the metastable levels of O 
VII and Ne IX, respectively, affecting the ratio just as a 
high density wou ld. We do not have s imultaneous UV 
observations, but iCollinge et al.l (|2001l) report fluxes of 



about 2 x 10~ 14 ergs cm -2 s _1 A -1 at both wavelengths. 
With oscillator strengths and A values from the CHI- 
ANTI database (jDere et al.ll2009[ ). we estimate that the 
UV excitation would be comparable to the radiative de- 
cay rates at distances of 1.8 x 10 15 cm (7200 R g ) and 
5.7 x 10 14 cm (2300 R g ) from the central source. 

While it is not possible to completely disentangle the 
density and distance, we can place some limits. The 
O VII ratio is consistent with the low density limit. It 
implies a density below 4 x 10 10 cm~ 3 and a distance 
from the central source above 1.8 x 10 15 cm (7200 R ff ) 
. The Ne IX ratio has larger uncertainties, but taken 
at face value it would imply a density between 6 x 10 
and 3 x 10 12 cm -3 or a distance from the center of about 
6 x 10 14 cm (2400 R g ). It seems somewhat unlikely that 
the Ne IX lines would be formed at a higher density and 
smaller distance than the O VII lines, and we conclude 
that one or both of the Ne IX lines is probably affected by 
the complex Fe absorption structure in that wavelength 
range. 

Additional sets of line detections may suggest a higher 
density region that is more clearly relevant to the out- 
flow. The previous emission lines were only assumed to 
originate in the same outflowing region as the absorption 
features. However, the detection of absorption features 
that also denote a similar density gives further support 
to this assumption. These absorption features are the 
Fe L-Shell density diagnostics of the Fe XXII absorption 
pair at 11.92 A and 11.72 A, and the Fe XXI absorp- 
tion pair at 12.285 and 12.38 A. The Fe XXII pair of 
lines has a much weaker detection than the triplet emis- 
sion lines, but we do find the 11.72 A line at a 99.997% 
confidence level, or 4er. The weaker of the two lines, at 
11.92 A, is detected at 2a confidence level, and we do 
expect this line to be present if the 11.72 A line is ob- 
served and collisional excitation dominates. The ratio of 
11.92 A to 11.72 A is 0.4 ± 0.2, and plac es the ratio in 
the lo w de nsity limit. Using t he work by iMauche et al.l 
((2001 and IMiller etHI (l2008h . we find an upper limit 
of7xl0 12 cm _3 . Furthermore in this density regime of 
n e < 7 x 10 12 cm ~ 3 , we would also expect the lines from 
Fe XXI. We do note that the Fe XXI line at 12.285 is de- 
tected at a 99.99% confidence level but not the Fe XXI 
12.38. See Figure[2] An upper limit to the normalization 
of the 12.38 is found to be -2.6 xlO -7 photons cm -2 s _1 . 
This places an upper limit to the ratio of 12.38/12.29 of 
0.1. Using Figure [71 this places an even stricter upper 
limit on this density region of n e < 6 x 10 11 cm~ 3 . 

2.4. The Photoionization Models 

2.4.1. Parameters 

Because the data are rich with information, we further 
characterized the X-ray absorption and emission features 
by implementing more self consistent and physical mod- 
els. To do so, we used two independent photoionization 
codes to create these models via XSTAR, version 2.2.0, 
and Cloudy, version 08.00. These photoionization codes 
are initially set up with a simple geometry including a 
spherical cloud, illuminated by a central source, i.e. an 
AGN. Beyond this, complexity is built into the mod- 
els by specifying the density, geometry and parameters 
of the illuminating source. This technique using multi- 
ple photoionization models allows for rigorous tests of 
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the components required by the data. However, a de- 
tailed comparison between the two packages is outside 
the scope of this paper. 

We begin with our average incident source, which has 
a ionizing luminosi ty of 10 42 ergs s" 1 , cons istent with 
iMaitra et all (|2Qlll ) and iBlustin et all (|2005f ). The inci- 
dent spectrum has a power-law shape in the X-ray regime 
with spectral index T = 2.0. Cloudy assumes a more 
detailed profile in other wavelength regimes, i.e. a black- 
body in the UV, while XSTAR only assumes a power-law 
from O.leV to MeV. However, due to lack of contempo- 
raneous data in other wavelength regimes, (e.g. in the 
UV) , a po wer-law descr i ption is adequate for this work. 
Moreover. IMaitra et al.l (|2011[ ) discuss evidence that the 
weak UV component in NGC 4051 may arise partially 
from a jet, rather than just the accretion disk. Evidence 
for jet prod uction in the centr al re gion of NGC 4051 is 
also given in King et al.1 (pOlTI ) and I Jones etall (pOlTI ). 

As the incident spectrum interacts with the absorbing 
material along our line of sight, the observed features 
can be produced. We derive the line widths (200 km/s) 
used in the simulations from the Gaussian fits to the 
absorption features. The covering fraction of the cloud 
is approximated at 75%, taken from the evidence of P- 
Cygni profiles , (i.e. Figure EB and aver age covering frac - 
tions given bv lBlustin et alJ(|2005D and I Reynold! (|1997t ). 
However, as will be shown below, the covering fraction 
only impacts the normalization of the emission lines. In 
addition, the inclination angle derived from the relativis- 
tic lines in the subsequent Section 12.4.31 suggests that 
the wind is neither a narrow wind along the polar axis 
nor a wind along the disk surface. Thus, a large covering 
factor is a reasonable assumption. A solar metallicity 
although not derived, was assumed for each model. 

Finally the density of each cloud was set according to 
the tentative density limits found from our line ratios in 
Section O A density of n e < 6 x 10 11 from Section OH 
used as a constraint on the initial set of models, which 
were created with densities of 10 10 and 10 11 cm -3 . The 
density of these regions is held constant, and we find 
that both density regions are required by the data. The 
10 11 cm -3 characterizes the more highly ionized regions, 
based on the conservative Fe XXII and Ne IX constraints, 
and the 10 10 cm -3 delineates the less ionized components 
at larger radii, based on the O VII triplet diagnostics. 

The data require three regions for both XSTAR and 
Cloudy models. The absorption models were included 
as multiplicative 3-dimensional tables varying the ion- 
ization parameter, column density and redshift for XS- 
TAR and inner radius, column density and redshift for 
Cloudy. The emission models were included as additive 
4-dimensional tables, which varied the above 3 parame- 
ters as well as a normalization parameter. We linked the 
column density and ionization parameter for the XSTAR 
absorption and emission components, and the column 
density and inner radius for the Cloudy absorption and 
emission components. Thus, the redshift and normaliza- 
tion were allowed to vary in the emission components. 
To compare between the parameters of the two models, 
we use the relation, 



where £ is the ionization parameter, Li on is the ionizing 
luminosity, n(r) is the electron density and r is the radius. 

2.4.2. Results 

Various combinations of the different density regions 
were tried, but inevitably two faster outflowing, highly 
ionized regions (log£ « 3.3 — 5.0) and a slower, less ion- 
ized region (log ^ « 1 — 2) were always present when glob- 
ally fitting the spectra for XSTAR and Cloudy. The best 
fit was in fact given by three absorption regions, two at 
10 11 and one at 10 10 cm~ 3 with a x 2 jv = 8177/7516 and 
X 2 /v = 8713/7516 for the XSTAR and Cloudy models, 
respectively (see Table 0] and [3]) . For completeness, we 
found a formally worse fit with three absorbing regions, 
at a density of 10 10 cm" 3 , with a x 2 jv = 8205/7516 for 
the XSTAR model. Assuming this is a change of two 
degrees of freedom, i.e. the density, and using an Ftest, 
we find that the XSTAR model with 10 11 cm -3 compo- 
nents is a 4cr improvement over models with 10 10 cm~ 3 
components. Repeating this with the Cloudy models, we 
found a marginally better fit with three 10 10 cm~ 3 at 
X 2 jv = 8709/7516. However, this is not statistically sig- 
nificant at only a confidence level of 1.2 a. We decided 
to include the Cloudy models with two 10 11 cm -3 and 
one 10 10 cm -3 density components for like comparison 
with the XSTAR model. (See Figures [p]] and [Tl]) . 

In addition to using the Ftest to justify our use of the 
higher density models, we also refit the spectra with and 
without the Fe XXII and Fe XXI lines in an effort to 
assess the extent to which these density sensitive lines 
were driving the fits. This was done only in the XSTAR 
models. By removing 0.3 keV regions around the Fe diag- 
nostic lines, and refitting the data, we find the ionization 
parameters of all three components decrease from log(£) 
= 4.5, 3.3 and 1.0 to log(£)= 3.6, 2.6 and 0.8. However, 
if we only include the regions around the Fe diagnostic 
regions in the spectral fit, the ionization parameters of 
the two highest regions stay roughly the same, while the 
lowest ionization component increases; log(£) = 4.5, 3.6 
and 2.4. This is demonstrating that the Fe L-Shell diag- 
nostic lines are in fact driving the spectral fits, for when 
they are included the ionization stays the roughly same. 
Conversely, when the lines are excluded, the ionization 
of all three components tends to decrease to lower ion- 
izations. It therefore seems a reasonable assumption to 
use the higher density models. We have repeated this 
process assuming all the components had a lower density 
of 10 10 cm -3 , and find the same trend in the ionization 
parameters but a statistically worse fit for each scenario. 

As a caveat, one should be careful interpreting these 
results, as the densities are upper limits. Reducing the 
density will tend to increase the radius at which these 
features are detected. 

Besides density, we also explored the effects of velocity 
broadening and covering fraction. In order to do so, we 
varied these parameters in the best-fit XSTAR model. 
Due to the resolution of the HETG spectrograph, we do 
not expect to resolve lines widths less than 150 km/s. 
However, larger velocity widths can be resolved. We 
fit the data with an XSTAR model that has a velocity 
broadening of 300 km/s and find a worse fit at x 2 l v = 
8275/7516. Finally, we also verified our choice of cover- 
ing fraction. Reducing the covering fraction to only 25% 
in the XSTAR model results in a worse fit of x 2 l v — 
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8188/7516. Although the 75% covering fraction model is 
a 3 a improvement over the 25% model, the change in 
parameters is only in the normalization of the emission 
lines. 

In general, the absorption regions specified by the fits 
are quite similar between the two packages. The Cloudy 
models are typically outflowing at a higher velocity than 
the equivalent XSTAR component, but in both the XS- 
TAR and Cloudy models, the emission features are con- 
sistent with being redshifted with respect to their ab- 
sorption counterparts. Furthermore, in both the XS- 
TAR and Cloudy models, the highly ionized region of 
log£ ~ 3.3 — 3.6, predominately characterizes the higher 
energy absorption features. This includes the lines es- 
pecially below 14A such as the He-like Si XII and Mg 
XI, the H-like Mg XII, and the Lyman Series of Ne X. 
Conversely, the lowest ionization region predominately 
characterizes the lower energy absorption features. This 
includes the He-like O VII lines and H-like O VIII lines 
above 16 A. This log£ « 1-2 also characterizes the O V 
feature at 22.4 A. Additionally. iCollinge et alj(|2001l) also 
find a C IV absorption feature at 1548.2 A and 1550.8 
A, which is formed at log£ ps 1.4. 

Of the absorbing regions, the most highly ionized com- 
ponent is the least constrained for it reaches the maxi- 
mum ionization parameter in both of the XSTAR and 
Cloudy grids. Thus, the £ is a lower limit in this model. 
However this component is required by an F-test and 
is required to characterize the highly ionized Fe L-Shell 
lines including the density dependent, Fe XXII. The two 
highest ionization parameters in both models, and there- 
fore smallest inner radii, place the absorbing regions 
within a few hundred gravitational radii of the SMBH 
assuming the density is constant with radius. At this 
radius one would expect line broadening to be approx- 
imately the orbital speed, resulting in a broadening of 
ps 4xl0 4 km s _1 . Clearly this is at least two orders of 
magnitude higher than the observed 190 ± 150 km s _1 
line widths derived from our Gaussian fits. However, 
due to the geometrical effects of encoding orbital motion 
into the absorbing region because of a small illuminating 
source assumed to be at <10 R g and a distance of 500 
Kg, the resulting broadening would only encode 2% of 
the orbital motion, i.e. « 300 kms -1 . This is consistent 
with the ltr error of our measurements. Alternatively, if 
the density of these regions is not constant but scales as 
r -2 , as thought to be typical of winds (e.g.. [Miller et all 
2008), or even r , then the inner radius need not be as 
close to the SMBH and can also remedy this discrepancy. 

Where the XSTAR and Cloudy models differ seems to 
be predominantly in the range of 10 - 16 A, and they 
differ mainly in the log£ ss 1 — 2 components. As seen 
in Figures [TU] and [TTJ one can readily notice the differ- 
ences in the absorption features especially around 14A, 
which is likely the Fe L- Shell charge states. The absorp- 
tion discrepancy between the two models is likely to be 
explained by the inclusion of the Fe L-shell transitions 
of Fe XVII - Fe XXII. The Fe L-Shell transition prob- 
abilities, ionization rates and cross-sections are less well 
known than the H-like and He-like transitions of O or 
Ne. So it should not be surprising the two models differ. 
The exact differences between the models are outside the 
scope of this paper, but we note that in a general sense 



the models follow the same trends. 

2.4.3. Fluorescence Fe Ka and O VIII Lyman a 

Also interesting in these fits are the lack of predicted Fe 
Ka or (3 emission at 1.94 and 1.75 A respectively. Cloudy 
fits the weak Fe XXIV emission line at 1.85 A but not the 
fluorescence lines at either 1.77 or 1.95A, and XSTAR 
only fits the slightest of emission at the Fe Ka wave- 
length. Conversely, both these lines as well as a broad 
component at a minimum of a 99.5% confidence level are 
required by the phenomenological fits. This is likely due 
to the fact that these fluorescence lines arise from neutral 
or very low ionization states of iron, and would not be 
modeled by the high ionization parameters in our models. 
We argue that these lines come from a different region 
entirely because of their detection in the phenomenologi- 
cal model and lack thereof in the photoionization models 
as well as the ionized states. The narrow lines are the 
reflection Fe K a and j3 lines probably produced in the 
outer regions of the accretion disk. Conversely, the broad 
component is thought to originate within a few hundred 
gravitatio nal radii of the blac k hole. This is consistent 
with what lYoung et all ()2005f ) find while studying MCG 
-6-30-15. They also state that the curvature produced 
by photoionization models could in principle produce a 
redwing effect in the Fe Ka regions, but we do not see 
this effect. 

As a result, we decided to model the Fe Ka and /3 
lines using two narrow Gaussian components fit to the 
data at 6.41 and 6.97 keV in the rest frame of the galaxy 
to describe the narrow components as well as Kerrdisk 
to model the broad component. We also include the XS- 
TAR models. Kerrdisk is a ray tracing model that ac- 
counts for both relativistic beaming and light bending 
effects in the vicinity of the acc retion disk and black hole 
(jBrenneman fe Reynolds! [2006D . In addition, it also as- 
sumes material within the radius of marginal stability 
is very highly ionized and has a low optical depth, thus 
lowering its contributions to the line profile. Finally, it 
further excels over other models, for not only does it 
allow the inner and outer radii, inclination angle, and 
emissivity index of both the inner and outer disk to vary, 
but it allows the spin parameter to vary. 

In our model, we fixed several of the parameters of the 
Kerrdisk model including the energy of the line to 6.39 
keV in the rest frame of NGC 4051, the emissivity index 
of the disk to be 3, the inner radius to be 1 R ms , and 
the outer radius to be 400 R ms , where R ms is the radius 
of marginal stability. We find that the best fit values 
using Kerrdisk at a x 2 /v = 8085/7507 and the model 
is plotted in Figure [8] We find a normalization of 4.0 ± 
0.6 xlO -5 photons cm -2 s _1 , and a maximum spin of 
a=0. 998_o. 088 with an inclination angle of 31.6l^'y°. A 
maximall y spinning b l ack h ole is in agreement with the 
results of lOgle et al.l (12004T) and the inn er radius is in 
agreement with iSteenbrugge et al.l ((2009). However, the 
inclination we find is lower and inconsistent with the 48° 
that both report. 

In addition to the Fe Ka emission line, we also explored 
the poss ibility of a relativ istic O VIII Lyman a line foun d 
by both lQgle ■ et~all (|2004l ) and ISteenbrugge et al.l (j2009fl . 
Evidence of such a line would suggest an over- abundance 
of Oxygen, which may be similar to the over-abundance 



G 



of Ir on seen in Seyfert-1, 1H 0707-495 (jFabian et al.l 
2009). We include this line with Kerrdisk fixing the 
rest frame energy at 0.654 keV {x 2 /v = 8169/7513). This 
component was required at a confidence level of just over 
2a, at 97.3%. The resulting parameters of the best fit 
model include a spin of a=0.998 (unconstrained) with an 
inclination of 37.6°t 4 4 and a normalization of 1.4 

xlO photons cm~ 2 s . Again we find a maximally 
spinning black hole, although not constrained, consistent 
with the aforementioned Fe Ka line parameters. The in- 
clination angle ha s incre ased but is not c onsist ent with 
iSteenbrugge et al.l <|2009T ) and lOgle etaLl (|2004| 1. 

If we tie the spin parameter, inclination angle and inner 
radius of the O VIII Lyman a line to that of the Fe K 
a and then refit the data, we find a similar significant 
detection at just over 2a, 96.9% {^jv = 8080/7506) . 
However, the Fe Ka line does drive the spin to 0.998_o.os 
with an inclination angle of 31.7+q'7°. The normalization 

of the O VIII Lyman a line is 1.2 ± 0.5 xl0~ 4 photons 
cm~ 2 s , and the normalization for the Fe K a line is 
4.0 ±0.6 xlO -5 photons cm -2 s _1 . This model, along 
with the aforementioned models is plotted in Figure |H1 

3. DISCUSSION 

3.1. Physical Characteristics 

Finding three distinct ionizing regions illuminates dif- 
ferent attributes of the winds in NGC 4051. We are able 
to constrain where these components are launched, the 
size of the emitting regions, how ionized the regions are, 
and the total amount of material that is observed to leave 
the system. 

Table [S] lists the aforementioned ionization parameters 
and radii of each of the three regions det ermined via XS- 
TAR and Cloudy. iDennev et all §009) use reverbera- 
tion mapping techniques to place the Broad Line Region 
(BLR) at 1.87 ± 0.5 light days or 4.84 ± 1.4 xlO 15 cm. 
The lowest XSTAR ionization component is consistent 
with this distance, while all the other components are in- 
terior to the BLR. This proximity to the SMBH suggests 
that the accretion disk may play a role in the launching 
mechanisms. 

Also listed in Table [6] is the radial extent of each re- 
gion. This is estimated by using the simple approxima- 
tion A(r) = Nji/n. The size of each region is between 
2 xlO 9 and 2.2 x 10 10 cm, less than 1% of the radial 
distance to the source. 

Because each region has a relatively small radial ex- 
tent, the mas s outfl ow rate is low. Following the work of 
IBlustin et all (pool , the mass outflow rate is determined 
assum ing a modified spherical outflow. As IBlustin et al.l 
(2005) describe, the outflow is most likely in clouds and 
filaments, so caution should be taken when interpreting 
these results. The outflow rate is described as follows, 

M out =r 2 p(r)v(r)n w 1.23r 2 m p n(r)vtl (3) 

Here we assume solar abundances, p— 1.23 m p x n(r), 
where m p is the mass of a proton, and a covering fraction 
of fl= 0.75 x 4-7T. It should be noted, that the density is 
assumed to drop off as 1/r 2 , so n(r)=n(R)R 2 /r 2 . To re- 
late the average number density, n(r) to the measured 
number density, n(r), a volume filling factor, C v (r), 



needs to be included as so, 

n(R) =n(R)C v (R) 



(4) 



A volume filling factor C„ is derived for each case using 
the following form, 



£N„R 



(5) 



where £ is the ionization parameter, N# is the column 
density, R is the distance from the source to the absorb- 
ing region and Li on is the ionizing luminosity. These 
values are given in Table |H1 
Finally, combining equations [21 El and 0] one gets, 



M n 



l.23m p Li on C v (R)vn 



(6) 



The resulting mass outflow rates are given in Table [6j 



+9.6 



The total mass outflow predicted by XSTAR is 7.9I" 
10~ 6 Mgyr -1 . Cloudy models give a higher predicted 
outflow rate of 1.3 + S | x 10" 5 M P 



In both models, 



the less ionized and slower moving component dominates 
this result by almost an order of magnitude. This can be 
attributed to the large covering factor assumed to be 75% 
and the large radius at which it is observed. The outflow 
rates may be lower limits, if gas flows at smaller radii 
with a very high ionization that inhibits its detection 
using atomic features. 

The observed mass outflow rate is very low compared 
to the implied mass accretion rate onto the SMBH. The 
mass accretion rate is, M acc = L^ i/(c rf), where Lf, ; is 
the bolo metric luminosity of the NGC 4051(L boi =10 43 - 4 
ergs s -1 . IBlustin et ai1 [2"005). and 77 is the efficiency fac- 
tor, typically assumed to be 10%. With M acc = 0.0047 
M0j/r _1 , one finds that M l °^ 1 /M acc is approximately 
0.002 and 0.003 for XSTAR and Cloudy respectively. For 
reference, the Eddington accretion rate would be MEdd ~ 
0.037 MQ-yr^ 1 , a much higher rate than either mass ac- 
cretion rate or the outflow rate. The low mass outflow 
rate can be attributed t o the very small fil ling factors 
derived from equation [5l IBlustin et al.1 ()2005| ) also find a 
small mass outflow rate of 0.0008 MQyr" 1 , but it is still 
an order of magnitude larger than the mass outflow rates 
found in this paper. 

3.2. Physical Scenario 

In addition to the amount of material that is be- 
ing driven from the central regions, we can also com- 
ment on the physical nature of this system. In previous 
work, authors have claimed evidence for both shocks , 
and two phase regio ns (i.e.. iPounds &; Vaughanl [20111 : 
iKrongold et~aLll2007f ). Here we discuss and compare our 
work to previous analysis of the warm absorbing winds 
in NGC 4051. 

The three regions in both of the XSTAR and Cloudy 
models follow the trend where the most ionized regions 
are closest to the source and fastest moving. T his is 
in contrast to the work done by lNucita et al.l (|2010l ) who 
find only one absorption component, which is moderately 
ionized and not outflowing. The lack of blue shift may 
be result of the lower r esolu tion of the RGS spectra. In 
addition. iNucita eT al. (2010) also find two emission com- 
ponents, one of which is consistent with our log£ = 3.3. 
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On the other h and, lLobban et al.l (|2011f) and 
iPounds fc Vaughanl (|2011l ) find a similar trend to 
our work, where the outflowing velocity decreases 
with decreasing ion ization parameter. In particular, 
lLobban et all |2011) find four zones with outflowing ve- 
locities ranging from -180 to - 820 km s _1 and ionization 
parameters ranging from log£ = -0.86 to 2.97. They also 
find an additional highly ionized region, log£ = 4.1 which 
is outflowing very quickly at 5800 km s . We do not find 
the lowest outflowing nor the lowest ionization regions, 
but our data are consistent with all but the largest veloc- 
ities, i.e. 400 < v out < 710 km s" 1 . This is likely due to 
the fact that the the HETGS is less sensitive to the low 
ionization features. Because there is a trend between in- 
creasing velocity and ionization, we would not expect to 
detect these slow moving, low ionization features. Table 
[71 lists these paramete r s for comparison. 

IPounds fc Vaughanl (|201lD find five absorption regions 
with velocities ranging from +120 to -5880 km s^ 1 with 
ionization parameters spanning log£ = 0.32 to 2.97. 
Again, we do not find the lowest velocity component 
nor the two highest velocity components, which are out- 
flowing faster than 3800 km s" 1 (See Table [7J. Sev- 
eral of their high velocity components are outside our 
wavelength range, but a majority of these highly ion- 
ized, highly velocity components are within our range. 
W e note that the componen ts b and c seen Figure 2 
of IPounds fc Vaughanl (|2011l ) are also readily explained 
by an O VII line rather than a highly blue-shifted 
O VIII Lyman a line, as shown in our XSTAR and 
Cloudy models (see Figur es [TP1 and [TT]). Additionally, 
IPounds fc Vaughanl (|2011l ) note that the high velocity 
components they observe with the EPIC CCD cameras 
include Mg XII through Ar XVIII. We do not find these 
in our high resolution spectra. 

T he multiple regions d e scribe d in lLobban et al.l (|201 ID 
and IPounds fc Vaughanl (|2011f ) are both discussed as a 
continuous region that spans a variety of different ion- 
ization states in the vicinity of a shocked shell. The au- 
thors suggest that their multiple regions are consistent 
with a constant mass outflow rate, i.e. M ut oc Lio ™ v = 

constant. Figure 8 in IPounds fc Vaughanl (|2011l ) demon- 
strates the necessary linear dependence of the observed 
ionization parameter and velocity. However, we do not 
find a similar linear dependence between the velocity and 
ionization parameters. Our ionization parameters span 
four orders of magnitude while the velocity shifts only 

s pan one. 

iSteenbrugge et all (|20()9: also find a number of ion- 
ization regions in their a variability study of NGC 4051 
using Chandra LETGS. In their highest signal to noise 
spectrum, C, they note four regions with ionization pa- 
rameters ranging from log£ = 0.07 to 3.19. Ins t ead o f 
evidence for a shocked flow, ISteenbrugge et all (2009) 
state that the ionized regions must originate within a 
0.02 - 1 p c distance of the central engi ne, i.e. much 
closer than IPounds fc Vaughanl (|2011l) and lLobban et al.l 
(poll . This distance is more consistent with our results 
but st ill formally discrepant. Further. ISteenbrugge et al~l 
(2009) do note a lack of variation in three ionizing ab- 
sorbers, placing a lower limit on the distance to the 
source, we suggest that there is a different scenario that 
could account for the lack of variability. If there is a 



continuous, flat distribution of absorbing column den- 
sity as a function of ionization parameter, a change in 
the ionizing flux could shift the ionization parameters in 
such a way that there would be little change in the ab- 
sorption lines; likewise, one would not see the associated 
variability following a change in the continuum. This 
allows the material to lie closer than previously derived 
from reionization timescales. We note that one compo- 
nent does show slight variability, component 2, which 
shows a 2.3er drop in ionization parameter over the four 
epochs. This is di s crepan t with this scenario. However, 
ISteenbrugge et all (|2009T ) note that they freeze the col- 
umn density for the spectra with low signal to noise, 
which could influence the ionization parameters because 
the two parameters are not clearly independent of one 
another. This may influence the aforementioned vari- 
ability. 

In addition, ISteenbrugge et all (|2009l ) find a range of 
velocities for their absorbing regions, from -210 to -4670 
km s _1 . Again, we do not find that our velocities are in 
agreement with the highest outflows. This is the result 
that our highest ionizing component, log £ > 4.5 charac- 
terizes the Fe L-shell transitions, instead of the C VI and 
N VII Lyman a lines. The velocity as well as ionization 
parameters al so given in Table [TJ for easy comparison. 

Conversely, iKrongold et al.l (|2007D find a two phase 
medium which is not consist e nt wi th t he shocks de- 
scribe d by IPounds fc Vaughanl (|2011f ) and lLobban et a~ 
(1201 If) or the multiple regions given by Stcenbrugge et al. 
(2009). Instead, the authors conclude that the two ab- 
sorbing regions are in pressure equilibrium with a pres- 
sure on order of 10 12 K cm 3 . A constant pressure would 
conflict with a constant mass outflow rate. In our work, 
we constrain the temperature to be 5xl0 5 K in XSTAR, 
forcing the two highly ionized regions to be in pressure 
equilibrium, and the third region not in equilibrium. 
Cloudy finds the equilibrium temperature which varies 
slightly between the models resulting in temperatures of 
approximately 4xl0 4 , 3xl0 5 and 7xl0 5 K for the com- 
ponents of logn e = 10, 11 and 11 respectively. Because 
the Cloudy fit is poor, there is only suggestive evidence 
that the components are not in pressure equilibrium. We 
would expect the temperature to increase with decreas- 
i ng density for a constan t pressure. 

IKrongold et al.l i|2()()7 l also find a lower electron num- 
ber density than we do, by 3-4 orders of magnitude, using 
variability analysis of an approximately 50 ksec XMM- 
Newton RGS spectrum. They also find higher ioniza- 
tion parameters than we find, and their photoionization 
model fails to characterize particular absorption features, 
most notably the Fe XVII line at 15 A. (See TableQ This 
discrepancy may be the result of the Chandra HETG 
spectra used in our work, which affords us a better sen- 
sitivity. However, we would like to note a s imilar trend 
between our work and IKrongold et ail (|2007l ). in that we 
find at least two regions with ionization parameters that 
differ by more than two orders of magnitude. 

In general we can make these conclusions following a 
detailed comparison to the literature, 

1. NGC 4051 is highly variable and discrepancies be- 
tween different observations are to be expected. 
Moreover, the different resolutions and wavelength 
ranges of XMM and Chandra will emphasize dif- 



ferent features. 

2. We find a low ionization component between 
log£ = 1 — 3 that is blue-shifted at rate of 200- 
600 km s _1 consistent with other work on NGC 



20111 ) 



4051 (e .g..lKrongold et al.ll2007t iSteenbrugge et al 



20091: lLobban et all 120111 : iPounds fe Vau ghan 



3. We do not find c ompon e nts log£ < 1 , unlik e 
ISteenbrugge et all (120091) . lLobban et al.l (j2011l ). 
and IPounds fe Vaughanl (|201lD ~ 

4. We do not find the high vel ocity compo- 
nents , i.e. > 2000 km s -1 (iCollinge et all 
20011: ISteenbrugge et~aTI 120091 : l Lobban et a.l.l l2011b 
Pounds fe Vaughanl 120111 ) . These high velocity, 
highly ionized compon ents characterize C VI and 
N VII Lyman a lines (Ste enbrugge et al . 2009) as 
well t he O VIII Lyman a lines ((Pounds fe Vaughanl 
l20ll and the Fe XXV and Fe XXVI lines at 6.8 
and 7.1 keV respectively (jLobban et al.ll2011[ ). The 
velocity shifts that these authors find, although 
high, are inconsistent with one another (See Table 
[7|) . Our work shows the highest ionization parame- 
ters are required by the Fe L-shell charge states but 
not at such velocities. Further more, the high veloc- 
i ty co mponent c in Figure 2 of IPounds fe Vau ghan 
([20T1 . the O VIII Lyman a line can also be in- 
terpreted as the O VII 18.627 A absorption line. 
The ls2 - ls3p line at 18.6 A is expected to be 
1/4 as strong as t h e 21. 6 A Is2-ls2p line, but 
IPounds fe Vaughanl (|2011l ) Figures 2 & 3 show it 
to be approximately the same strength. We inter- 
pret this as partial filling of the 21.6 A absorption 
feature by the resonance line emission that must 
be present along with the observed forbidden and 
intercombination lines. However, this could also be 
due saturation of the 21.6 A line. 



5. ISteenbrugge et all (|2009| ) discuss lower limits to 
radii due to lack of variability in the ionization pa- 
rameters; we suggest that a continuous set of ion- 
ization parameters could also account for the lack 
of variability at a smaller radii. 



6. Follow ing I Ogle et"al"1 ([2004D and ISteenbrugge et al.l 
(2009), we can also report the possibility of an O 
VII Lyman a relativistic line. In addition, we find 
the data require a Fe Ka relativistic line. Taken lit- 
erally, the fits suggest a maximally spinning black 
hole, a=0. 998_o. 08 with an inclination angle to 
31.7^q'7 , when these parameters are linked be- 
tween these two relativistic lines. 

Our analysis shows characteristic traits in common 
with both the two phase medium, i.e., at least two com- 
ponents that span a wide range of ionization parame- 
ters, and the continuous range in ionization parameters 
and velocities. However, neither of these scenarios are 
entirely consistent with the data treated in this paper. 
As a result, we believe that we may be seeing at least 
two regions separated by a large distance. The inner- 
most absorbing region is within a few hundred R 9 and is 
made up of at least two highly ionized components, i.e. 



log ^ > 3.3. which may be in pressure equilibrium. In 
fact, there may even be a range of ionization parameters 
in this region. The other component which is further out 
is also outflowing, but may not be related to the inner 
part of the accretion flow. The outer component is not in 
pressure equilibrium and could be slowing down due to 
inter actions with th e host galaxy. It is interesting to note 
that iDetmers et~al"1 (|201lD describe in Seyfert Mrk 509, 
a two-component absorbing region with high ionization 
parameters, and a te ntative third lower ioni zation region. 
Similar to this work, IDetmers et al.l (|2011l ) describe two 
discrete components which show a range in ionization 
states within each component but are not continuous be- 
tween the two. In addition, the two faster outflowing 
components at -319 and -77 km s" 1 found in M rk 509 
are also similar to our work (jDetmers et al.ll2011l ). 

3.3. Launching Mechanisms 

Regardless of the exact geometry, it is quite intriguing 
that both of the photoionization codes place a portion 
of the absorbing material very close to the central en- 
gine, suggesting that it is launched in the vicinity of the 
accretion disk. At such a small radii, thermal pressure 
is ruled out as a mechanism for driving the winds. The 
escape temperature at a few hundred gravitational radii 
is greater than 10 11 K, while the gas temperature in our 
models is only on order of 5 x 10 5 K making escape due 
to thermal motions quite difficult. Consequently, either 
radiation pressure or magnetic fields are needed to drive 
the winds. To be driven by radiation pressure, the winds 
need to have a strong source of UV radiation as well as 
a large cross section to interact with the UV photons. 
However, each of these absorbing regions has a ioniza- 
tion above log£ > 3.3, which lower s the radiation pres- 
sure multiplier to very close to 1 (|Kallman fe McCray 
1982). Therefore, the gas does not gain any significant 
additional acceleration fro m radiation pressu re at these 
high ionization states (e.g.. iProga et al.ll2000D . 

This leaves the most probable launching mechanism 
to be mag netic fields. One particu lar scenario is de- 
scribed by iBla ndford fe Pavnd (|1982l ) who suggest that 
the winds can be launched from magnetic field lines that 
are tied to the accretion disk. The ionized gas flows 
along the magnetic field lines, and due to the rotation 
of the disk, centrifugal forces launch the ionized mate- 
rial. In this picture, the wind would carry away most 
of the angular momentum and much of the energy of 
the accreting gas. Therefore, the standard estimates 
for mass accretion rate based on luminosities would be 
an underestimate. Another way to launch the winds 
via magnetic fields is through magnetic pressure gener- 
ated through the magneto rotational instabilities (MRI) . 
iProga fe Begelmanl (|2003l ) predict this type of wind to 
generally have high density and low velocity, which is in 
accordance with our observations. One of the main dif- 
ferences in these two models is that the centrifugal winds 
are launched vertically where as the magnetic pressure 
winds are launched torodially. Unfortunately, we do not 
have enough sensitivity to measure any rotational veloc- 
ity encoded into the absorption features to distinguish 
b etween t he tw o. 

iBeharl (|2009f ) also finds evidence for magnetically 
driven Seyfert outflows, which lends further support to 
our observations. He performs a study using the ab- 



9 



sorption measure distribution (AMD) as compared to the The authors would like to thank the anonymous referee 

ionization parameter, log£, of 5 Seyferts. In general, he for their invaluable comments and improvements to this 

finds a very shallow relation between AMD and log£. paper. 

This can be interpreted as consistent with non-spherical 

MHD outflows, assuming the variations in log£ are not 

due to a turbulent interstellar medium. The resulting 

density profile scales inversely with radius, n oc r _1 . 



4. CONCLUSION 

We have used two separate modeling techniques, a 
phenomenological model and photoionization models, to 
quantify the absorption and emission features seen in the 
Chandra HETG spectra of NGC 4051. Amongst the re- 
sults are three sets of line detections that afford us the 
ability to diagnose the density of the material. These 
are the He-like Ne IX and O VII triplets as well as the 
tentative Fe L-shell doublet, Fe XXII. We also constrain 
an upper limit to the Fe XXI ratio, another Fe L-shell 
density diagnostic. In particular, the fact that the Fe 
XXI and Fe XXII ratios are absorption features means 
that they are particularly effective probes of the density 
of the absorbing regions in so much that it is created 
within the absorbing region itself, whereas the emission 
features may come from another nearby region. 

Following this phenomenological study, we imple- 
mented more self-consistent photoionization models via 
XSTAR and Cloudy. We used electron densities of 
n e = 10 10 and n e = 10 11 cm~ 3 based on the gas density 
diagnostics, and the data then require three different ab- 
sorption regions. The two highest ionized regions with 
log £ > 3.3, are at the higher density, while the lower den- 
sity region is less ionized, with log£ ~ 1 — 2. All three 
components may originate just at or within the BLR. 
The highly ionized regions likely originate close to the 
black hole, i.e. less than 300i? g (assuming no radial drop- 
off in density) and have a large velocity as compared to 
the lower ionization region that has a lower velocity and 
resides farther out, i.e. greater than 3300i? s . Velocity 
broadening of the inner-most components would nomi- 
nally show larger absorption line widths. However, in 
Section 12.4.21 we demonstrate that a small central source 
(<10 Kg Char tas et al.l I2009T ) may only impart a small 
portion of the rotational velocity in the absorption com- 
ponents. 

Furthermore, we demonstrate that the mass outflow 
rate is small compared to the overall mass accretion 
onto the black hole. It is possible that magnetic 
fields are partly responsible for driving these winds, at 
least in the regions closest to the black hole. This is 
very similar to th e galactic bl ack holes, GRO 165 5-40 
(iMiller et al.1120011. H17 43-322 (j Miller et al.ll2006b[) . 4U 
1630 -472 dKubota et al1l2007D th e dwarf nova, OY Cari- 
nae (iMauche fe Ravmond[ l2000^ and the Seyfert, NGC 
4151 (jKraemer et al.ll2005t l. in which magnetic fields are 
also found to potentially be important in driving out- 
flows. 

Finally, we note the data require a relativistic Fe Ka 
line, which nominally implies a maximal spin 0.998_o.O8 
with an inclination angle of 31.7^07°- At a 2 a level, 
the data require a rela tivist ic O VIII Lyma n a line. 
iSteenbrugge et all (|2009D and lQgle etahl (|2004ft also find 
evidence for a relativistic O VIII Lyman a line. 
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TABLE 1 
Broad-Band Parameters 



r 


Normalization 
(10~ 3 photons cm" 2 s" 1 kcV" 1 ) 


Constant 


Diskbb 
(kcV) 


Normalizatoin 


1.94 ± 0.01 


4.24 ± 0.03 


1.062 ± 0.005 


0.118 ± 0.005 


6610 ±1™ 



Note. — This table gives the parameters for the broad-band continuum of both the HEG and MEG spectra. Because of a soft excess 
component in the MEG spectra, a disk blackbody component was included. The disk blackbody component should not be regarded as 
physical. The multiplicative constant is applied to the MEG spectrum to correct for the slight discrepancy between the flux in the two 
spectrographs. Galactic absorption of 1.15 xlO 20 atoms cm -2 was also included. X 2 /" = 9067/7531 



TABLE 2 
Emission Line Parameters 



line center 


a 


Normalization 


EW 


Ac 


velocity 




Ftest 




ID 


(A) 


(mA) 


(10~ 5 ph cm" 2 s" 1 ) 


(mA) 


(A) 


(km s _1 ) 










22.135 ± 0.003 


12.8 tii 


8.73 


132 ±21 


22.101 


-240 t1° 


5 


.97x10" 


-08 


O VII 


91 SSI +0-026 
Zl.oor _Q.027 


11.6 t^ 3 


2 54 + 1 - 30 
' ! ■ J, -0.88 


37.8 til;? 


21.807 


Q7n +360 
_370 


5. 


.37x10" 


-02 


O VII 


19.013 ±°;°gg 


12.4 ±fl 


9 t-n +0.79 
z ' oz -0.92 


37.6 ±ft» 


18.967 


IS +" 
18 -40 


2 


.10x10" 


-04 


O VIII 


18.668 


0.00 + 4 -' 17 


1.07 ± 0.31 


56.8 ± 16.6 


18.627 


40 + 31 
-100 +«o 


5 


.32x10" 


-03 


O VII 


16.038 ±°;g°t 


8.89 ttit 


1 oq +0.59 
L -° u -0.81 


48-0 till 


16.006 


3 


.84x10" 


-02 


O VIII 


13 716 + 003 
lO. 1 lO _Q.004 


8-34 tlH 


0.83 1°; 23 


17.9 +*■§ 


13.699 


-330 +™ 


4. 


.76x10" 


-05 


Ne IX 


13 579 +0- 005 
lo.oiy _o.oo6 


13.48 ±g-gi 


75 + - 18 
u - ' 3 -0.19 


16.3 +li 


13.553 


117 + 116 
1J -' -123 


1 


.15x10" 


-03 


Ne IX 


19 icr +0.005 
IZ.IOD _ .004 


64 + 8 ' 48 
U - D4 -0.64 


41 + 01 
u - 41 -0.13 


10.1 t°j 


12.134 


179 +113 
' Llz -92 


4. 


.04x10" 


-03 


Ne X 


S AR7 +0-001 
8-467 _o. 002 


07 + 3 - 10 
u - u ' -0.07 


0.21 ± 0.05 


5.7 ± 1.5 






1 


.80x10" 


-03 


(Mg XII ?) 


A 9SS +0-004 
4.28S _ .001 


0.00 + 5 - 87 


0.23 ± 0.07 


6.6 ± 1.9 


4.299 


1 /ico +290 

-1460 I 40 


2. 


.49x10" 


-03 


S XV 


1 OA1 +0.002 


1.48 t 3 ; 3 ° 


54 + - 17 

"J-04 _q 12 


10.8 t 3 ;| 


1.940 


-620 +23 ° 


4. 


.51x10" 


-04 


Fe Ka (narrow) 


1 Q97 +0-011 
l.yzi _o.oi2 


62.9 +{ti 


9 01 +0.43 
z - 0± -0.42 


65.9 +IH 


1.940 


-2720 tllfo 


1 


.22x10" 


-06 


Fe Ka (broad) 


1.781 ± 0.004 


10.00 ±tf 8 


n RA +0.16 

0.64 _ 17 


17.9 til 


1.750 


4530 t™ 


5 


.23x10" 


-03 


Fe K/3 



Note. — The above tabic lists the properties of the emission lines measured in HETG coadded spectrum of NGC 4051. A range of ionization 
and velocity shifts are measured. The first column lists the initial wavelength. The second and third list the width as well as the normalization for 
the corresponding line. Those lines with a a— are unresolved. The Equivalent width is given in the fourth column. Ao is the rest wavelength in 
NGC 4051 rest frame, which is at a z— 0.00234. The velocity is calculated using these rest wavelengths in column 6. The seventh column lists the 
significance of detection as 1 minus the confidence level. The final column lists particular ion that is being detected. All the line energies are fit to 
the data, which explains why the energy is different for the narrow and broad Fe Ka lines. 
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TABLE 3 
Absorption Line Parameters 



line center 




Normalization 


EW 


Ao 

(A) 


velocity 


Ftest 


N z 


ID 


(A) 


(raA) 


(l(r 5 ph cm- 2 s- 1 ) 


(mA) 


(km s- 1 ) 




(10 16 cm- 2 ) 





91 R17 +0 006 



+0.001 

-0.004 

-0.0012 
+0.003 
-0.008 



18.974 
18.650 
17.790 
16.020 _ 0074 
15.266 ± 0.004 
15.017 ± 0.002 
14.220 ± 0.006 
14 042 + 010 

_ .013 

13.841 ± 0.002 

13 658 +0- 003 
lo.ODo _ 005 

13.515 ± 0.002 
13.450 ± 0.002 
12.832 ± 0.003 

19 981 +0 003 
iz.zoi _ 004 

ir, lor, +0.012 

li,IJi -0.002 
11 935 +0 045 
i±.»ou _Q.0002 

11.771 ± 0.002 

+0.001 
-0.005 



11.655 



11.551 ± 0.003 
11 415 + 005 

ir.iro -0.0001 

11.001 ± 0.004 

m 9"37 +0.003 
lU.zof _ .004 
q SOS +0.002 
».ooo _Q.0003 

9 174 +0 003 

y - ±,yi -0.004 

8.421 ± 0.001 

6 650 +0 0° 24 
O.DOU _ 0018 



13.2 
15.7 
33.8 
1.29 
14.3 
5.38 
4.98 
28.2 
23.4 
8.41 
14.3 
8.18 
14.6 
9.78 
3.91 
9.12 
0.00 
2.90 
0.16 
6.11 
0.00 
8.03 
5.41 
0.00 
8.19 
2.44 
4.32 



+9.6 

-13.2 

+4.8 

-6.0 

+5.7 

-5.2 

+ 12.86 

-1.29 

+5.5 

-7.1 

+4.79 

-5.38 

+3.63 

-4.98 

+6.7 

-5.8 

+ 10.4 

-6.1 

+4.51 

-3.57 

+8.3 

-5.3 

+3.19 

-2.91 

+2.1 

-2.5 

+3.34 

-2.99 

+9.75 

-3.91 

+6.05 

-2.72 

+6.86 

+4.29 
-2.90 
+4.44 
-0.16 
+3.55 
-6.11 
+4.50 

+4.15 
-4.35 
+4.84 
-5.41 
+3.03 

+4.65 
-3.33 
+ 1.84 
-2.44 
+2.84 
-2.60 



-2.71 
-3.46 
-4.89 
-0.94 
-1.98 
-0.46 
-0.76 
-1.55 
-0.83 
-0.76 
-0.56 
-0.73 
-1.23 
-0.66 
-0.35 
-0.87 
-0.13 
-0.32 
-0.19 
-0.33 
-0.18 
-0.25 
-0.25 
-0.18 
-0.28 
-0.26 
-0.23 



+0.52 
-0.66 
+0.77 
-0.61 
+0.61 
-0.64 
+0.17 
-0.58 
+0.85 
-0.52 
+ 0.14 
-0.15 
+ 0.14 
-0.16 
+0.29 
-0.32 
+0.25 
-0.29 
+ 0.15 
-0.18 
+0.17 
-0.19 
+ 0.13 
-0.14 
+0.15 
-0.16 
+ 0.11 
-0.12 
+0.09 
-0.18 
+ 0.14 
-0.48 

± 0.07 

+0.07 
-0.08 

± 0.06 
± 0.08 
± 0.06 
± 0.08 

+0.07 
-0.08 

± 0.05 

+0.08 

-0.01 

± 0.05 

+0.06 
-0.07 



+7.6 
-9.7 
+ 10.7 



-39.9 
-48.2 
-70.8 
-14.8 
-34.5 
-9.0 ; 
-15.3 
-32.7 
-17.9 
-16.1 
-12.1 
-16.0 
-26.9 
-15.1 
-8.3 . 
-20.9 
-3.1 ± 1.7 



— 8.0 
+8.9 
-9.3 
+2.6 
-9.1 
+ 14.8 
-9.1 
-2.8 
■3.0 
+2.9 
-3.3 
+6.2 
-6.7 
+5.3 
-6.3 
+3.1 
-3.8 
+3.7 
-4.2 
+2.8 
-3.1 
+3.3 
-3.5 
+2.6 
-2.8 
+2.1 
-4.3 
+3.4 
-11.4 



-1.7 
-2.1 

-4.7 ± 1.6 
-8.4 ± 2.1 
-4.6 ± 1.5 

H.9 
-2.0 
hi. 6 
-1.9 

-4.8 ± 1.3 

+2.3 
-0.3 

-7.2 ± 1.4 

+ 1.6 
-2.0 



.1 



-6.3 
-6.4 



-7.7 



-6.5 



21.602 
18.967 
18.627 
17.768 
16.006 
15.261 
15.014 
14.203 

13.844 
13.670 
13.518 
13.447 
12.847 
12.285 
12.134 
11.920 
11.770 

11.547 
11.427 
11.001 
10.238 

9.169 
8.419 
6.648 



-497 
-584 
-331 
-332 
-438 
-592 



+79 
-111 
+ 10 
-63 

-19 
+43 
-126 



-138 
+77 
-76 

-641 ± 45 

hl32 
-137 



-334 



751 + 52 
-<o± _ 4 g 

Q71 +58 
" M,± -119 

-769 ±H 
-629 ±52 
-1046 ±H 
"779 ±f 06 
-737 ±lf 
-320 ±l 100 
-687 ±f 6 



-602 ± 74 
-1014 ±i 31 

-691 + 120 
DM± -121 
700 +85 

-(6S _ 103 



-529 +"° 
0ZJ -118 

-635 ±%l 

-612 + 109 
Diz -79 



3.04xl0-° 7 

8.29 XIO- 35 

8.37X10- 28 

6.11xl0-° 8 

3.94X10- 17 

4.17X10" 03 

1. 27x10-" 

2.06X10-" 8 

4.79X10" 03 

1.53X10" 09 

4.47X10" 03 

8.30x10-" 

6.48X10- 22 

2.33xl0- 1() 

1.07xl0-° 4 

7.62 xlO- 25 

4.6xl0-° 2 

3.54xl0-° 5 

4.40xl0-° 2 

l.OlxlO" 04 

3.48xl0-° 2 

1.13X10" 02 

5.60xl0-° 4 

3.53X10" 03 

6.97xl0-° 4 

2.23X10- 13 

3.34X10" 06 



1.2 ± 0.2 
4.8 ± 0.8 
14 ± 2 

8-5 til 
51 ± 8 
0.6 ± 0.1 
0.3 ± 0.1 
1.7 ± 0.3 

16.2 ± 2.6 
1.4 ± 0.2 
1.1 ± 0.2 

2.1 ± 0.3 

2.2 ± 0.4 
0.6 ± 0.1 
5.1 ± 0.8 
0.4 ± 0.1 
0.9 ± 0.1 



4 


2 


± 


0.7 





2 


± 


0.1 


9 


3 


± 


1.5 



11.5 ± 1.8 

1.2 ± 0.2 
3.7 ± 0.6 
1.9 ± 0.3 



O VII 
O VIII 
O VII 
O VII 
O VIII 
Fc XVII 
Fc XVII 
Fc XVIII 
Unidentified 
Fc XIX 
Fe XIX 
Fe XIX 
Nc IX 
Fc XX 
Fc XXI 
Nc X 
Fc XXII 
Fc XXII 
Unidentified 
Ne IX 
Fc XXII 
Nc IX 
Nc X 
(Fc XXI?) 
Mg XI 
Mg XII 
Si XIII 



Note. — This tabic lists the absorption lines detected in our phcnomcnological modeling in the HETG spectrum of NGC 4051. All of the features 
appear to be blue shifted with respect to the host galaxy, indicating a bulk outflow from the galaxy. The first column lists the initial wavelength. 
The second and third list the width as well as the normalization for the corresponding line. A negative normalization indicates an absorption 
feature. Those lines with a cr— are unresolved. The Equivalent width is given in the fourth column. Ao is the rest wavelength in NGC 4051 rest 
frame, which is at a z— 0.00234. The velocity is calculated using these rest wavelengths in column 6. The seventh column lists the significance of 
detection as 1 minus the confidence level. The eighth column lists the column densities derived using the relation W\ — (ire 2 /m e c 2 )NjX 2 fij , where 
W\ is the equivalent width, Nj is the column density, A is the rest wavelength, and fij is the oscillator strength. The final column lists particular 
ion that is being detected. All the line energies are fit to the data. 



TABLE 4 
XSTAR Model 



I Absorption | | Emission | 

Density Column Density log(£) velocity Normalization velocity A% 2 Ftcst 



(cm" 3 ) (10 20 cm- 2 ) (ergs cm s" 1 ) 


(km s- 1 ) 


(10- 3 ) (kms- 1 ) 




10" 8.10l^ 2 4.50_ 006 
10" 10.10I°,j 3 3.28±°'°| 
10 1() 1 -97i u '.25 l-00ig;?g 


-esol 2 " 

-630±i° 
-400l 38 ° 


19.53l 5 5 ; 22 8 -50±™ 

i.oiiS:S looii™ 

± - uo -0.16 OUU_ 40 


106 6.3xl0- 21 
1076 1.8xl0- 193 
248 26.5xl0 -48 


Power-law Normalization Constant 
(r) (10- 3 ph 

cm -2 s _1 kcV -1 ) 


Diskbb 

(kcV) 


Normalization Phabs 

((%7^) 2 -e) xlO 20 

(atoms cm - ) 




1.95 ± 0.01 4.32tj;;^ 1.059±£$2 


n , 00+0.002 

u.ioo_q 003 


3.78±£$ x 10 3 1.15 





Note. — The above table lists the parameters of the three XSTAR model components. The errors reported here are la. x 2 /v 8177/7516—1.09. 
The actual model is as follows: phabs*(mtablcl*mtablc2*mtablc3(powcr-law+diskbb)+atablcl + atablc2 -+- atablc3) :f! constant) The multiplicative 
constant is applied to the MEG spectrum to correct for the slight discrepancy between the flux in the two spectrographs. * Denotes the model 
reached a maximum in the ionization parameter limits 
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TABLE 5 
Cloudy Model 



| Absorption | | Emission | 

Density Column Density Radius velocity Normalization velocity A)? Ftcst 



(cm" 3 ) (10 2() cm~ 2 ) (10 13 cm) 


(km s 1 ) 


(10~ 21 ) (km s- 1 ) 




10 11 21.8+4 S 1.00+ 04 
10" 7.68t u; 8 75 5.09+_ il 
10 1() 6.42l°-j5 82.91H 


-1090±iJg 
-1000+™ 
-700± 30 , 


o.oo7± 0.002 -nolSo 

0.015± 0.006 220lgg° 
9 47 +0.80 9qn+ 60 

z - 4 '-0.7() ~ zyu -180 


23 1.4 x 10~ 3 
642 7.6 x 10~ 96 
484 1.9 x 10~ 54 


Power-law Normalization Constant 
(r) (10~ 3 photons 
cur 2 s _1 kcV -1 ) 


Diskbb 
(keV) 


Normalization Phabs 

(atoms cm -2 ) 




1.97 ± 0.01 4.42 +om lm ± °- 01 


n 197 -fU.UUa 
Uil ^' -0.003 


4.81 lo 57 X 10 3 1.15 





Note. — This table lists the Cloudy model parameters of the three components required by the data. The errors reported here arc la. x 2 / v 
—8713/7516—1.16. The actual model is as follows: phabs*(mtablcl*mtablc2*mtablc3(powcr-law+diskbb)+atablcl + atablc2 + atablc3)*constant) . 
The multiplicative constant is applied to the MEG spectrum to correct for the slight discrepancy between the flux in the two spectrographs. 



TABLE 6 
Photoionization Model 



XSTAR 



Cloudy 

2.16 ± 0.05 
82.9 ttl 
3310 t\ 9 7 
n fi4 +0 07 
U - D4 -0.08 

0.77 ± 0.13 

11 5 + 7 ' 1 
L1 -° -8.6 



XSTAR 



Cloudy 



XSTAR 



Cloudy 



Density (cm 3 ) 
log£ (ergs cm -1 ) 
Radius (10 13 cm) 
Radius (R 9 ) 
A(r) (10 10 cm) 

c v (10- 4 ) 

Mont (10- 6 M Q yr- 



10 1 



10 1 



4-50 '_ 0M 


5.00 


* 

-0.03 


3.28 


1.78 +°- 12 


1.00 


+0.04 
* 


7.25 


71 + 4 


40 + 


1 


290 


U-si _ . 10 


2.18 


+0.44 
-0.35 


1.01 


4-56 S fil 


21.8 




1.39 


0.30 + 8 ° 9 


0.73 


+3.28 
* 


1.42 



+0.05 
-0.03 
+0.29 
-0.38 

+ 11 

-15 
+0.06 
-0.09 
+0.17 

°+i 8 35 
-0.22 



3.59 ± 0.05 

5 09 + 29 
O.uy _ 26 

203 li, 1 , 
0.77 ± 0.08 

1-51 + °22 25 

1.20 if; 39 



1 00 + - 06 

iuu -0.15 
017 +58 
-20 

12700 till 

0.20 ± 0.03 

06 + 01 
+W 
-7.4 



10 10 



7.7 " 



Note. — This table compares the XSTAR and Cloudy model components. The models give similar values, and consistent mass outflow rates. 
We note that in general the lowest ionization component carries the most mass. We have assumed 1.7 ± 0.5 x 10 6 Mq as the mass of the SMBH. 
The symbol, *, denotes both upper and lower limits in the models. 
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TABLE 7 

Comparing Models 



Instrument 



Our XSTAR 
Model 
(HETG) 



Collingc 
ot al. (2001) 
(HETG) 



Krongold 
ot al. (2007) 

(RGS ) 



Stccnbruggc 
ct al. (2009) 
(LETGS) ** 



Nucita 
ct al. (2010) 
(RGS) 



Lobban 

et al. (2011) 

/"urn -1 /" 1 \ 
(Hhlb ) 


Pounds 
ot al. (2011) 

(Kbb j 


-0.86 ± 0.30 


0.32 ± 24 


-180 ± 100 


+120 ± 60 


_ 

< 23.4 


_ 
- 


20.48 ± 0.13 


20.0 ± 0.1 


n fin _|_ n Qfi 


i AO _|_ n on 


-220 ± 60 


-530 ± 90 


< 11 


~ 6.7 


15.8 - 22.6 


~ 18 


20.31 ± 0.13 


20.04 ± 0.26 


1.96 ± 0.17 


2.52 ± 0.07 


-820 ± 30 


-3850 ± 60 




~ 5.7 


< 20.3 


- 17.8 


20.70 ± 0.16 


21.26 ± 0.10 


2.16 ± 0.21 


2.77 ± 0.15 


-550 ± 60 


-400 ± 50 


< 20.4 




21.04 ± 0.27 


21.46 ± 0.29 


2.97 ± 0.10 


2.97 ± 0.06 


-710 ± 40 


-5880 ± 60 




~ 5.7 


^ 18-7 


- 17.8 


21.44 ± 0.12 


22.15 ± 0.07 


4.1 ± 0.2 




-5800 ± 1200 





log^o ( er g s cm s 1 ) 
v (km s -1 ) 
logno (cm -3 ) 
logro (cm) 
log N H ,o (cm" 2 ) 

log £i (ergs cm s ) 
vi (km s -1 ) 
log ni (cm -3 ) 
logri (cm) 
log Nh,i (cm -1 ) 

log 6a (ergs cm s -1 ) 
V2a (km s -1 ) 
logri2a (cm -3 ) 
logr2a (cm) 
log N H>2a (cm -2 ) 

log? 2i ,(ergs cm s -1 ) 
v 2b (km s -1 ) 
\ogn 2b (cm -3 ) 
logr 2i) (cm) 
log N H ,2b (cm -2 ) 

log £3 (ergs cm s -1 ) 
V3 (km s -1 ) 
logri3 (cm -3 ) 
logr3 (cm) 
log N Hj3 (cm -2 ) 

log £4 (ergs cm s -1 ) 
V4 (km s -1 ) 
logri4 (cm -3 ) 
logr4 (cm) 
log N Hj4 (cm -2 ) 



1.00 
-400 



+0.06 
-0.15 
+380 
-270 

10 

15.5 ± 0.1 
20.3 ± 0.1 



o no +0.05 

03 

-620 1™ 
11 

13.86 ± 0.02 
21.00 



+0.03 
-0.04 



^4.5 
-680 1 2 ° 
11 

13.25 + 03 

on 91 + - 05 
zu.ai _ .o6 



-600 ± 130 
< 10.6 

> 21.3 



-2340 ± 130 
< 10.6 

19 - 20.7 



2.18 ± 0.02 * 
-492 ± 97 
> 7.9 
< 14.9 
20.73 ± 0.17 



4.42 ± 0.07 * 
-537 ± 130 
6.8 - 7.3 
15.1 - 15.4 

21.42 ± 0.12 



0.07 ± 0.13 
-210 ± 70 

< 6.5 

> 17.95 
20.08 ± 0.18 

0.87 ± 0.09 
-200 ± 30 
> 6.6 
< 17.48 
20.46 ± 0.08 

2.32 ± 0.13 
-580 ± 50 

< 6.7 

> 16.70 
20.90 ± 0.3 



3.19 ± 0.09 
-4670 ± 150 
< 5.4 
> 16.85 
22.30 ± 0.3 



2.15 ± 0.02 

^ 7 
£ 16.8 
23.36 ± 0.01 



15.1 - 16.1 
22.92 ± 0.12 



Note. — This table shows the comparison between other absorption X-ray components found in the literature and the work presented in this 



paper. ♦Krongold ct al. (2007) does not use £ but instead uses the ionization parameter of U 



where Q(H) is the number of ionizing 



photons, r is the distance from the ionizing source, n e is the electron number density and c is the speed of light. To convert between the two, we 
assumed an ionizing luminosity of 10 42 . ** Used Spectrum C for this comparison. 
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Fig. 1. — This figure contains the ratio of the data to the fiducial model with x 2 / u 9607/7516. This model only includes a powerlaw 
and disk blackbody component as well as an effective H column density and a multiplicative constant, which corrects for the normalization 
between the HEG and MEG spectra. Clear absorption and emission features can be seen over this simple model. The wavelength plotted 
is in the observed reference frame. 
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Fig. 2. — These Ne X and Fe XXI lines are modeled with Gaussian 
profiles. Here one can see that the absorption is blue shifted with 
respect to the emission feature seen in the H-Like Lyman a line of 
Ne X. The P-Cygni-like profile of the H-like line could be real, but 
may also be a superposition of distinctive wind components. The 
line markers are in the rest frame of NGC 4051. The wavelength 
plotted is in the observed reference frame. 
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Fig. 3. — These OVII and OVIII lines are modeled with Gaussian 
profiles. Here one can see that the absorption is blue shifted with 
respect to the emission features seen in the H-Like Lyman series 
of O VIII. The P-Cygni-like profile of the H-like line could be real, 
but may also be a superposition of distinctive wind components. 
The line markers are in the rest frame of NGC 4051. The wave- 
length plotted is in the observed reference frame, and the data was 
rebinned with a minimum of 30 counts per bin for clarity. 
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Fig. 4. — This figure shows the OVII lines characterized with the 
phcnomcnological model for both absorption and emission. This 
He-like triplet gives a density estimate of n e as 10 10 cm -3 . The 
wavelength plotted is in the observed reference frame, while the 
line identifications are in the rest frame of the host galaxy. 
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Fig. 5. — This figure shows the Ne IX line complex that is mod- 
eled with Gaussian lines for both the absorption and emission. This 
He-like triplet gives a density constraint of n e <10 12 cm -3 . Line 
identifications are in the rest frame of NGC 4051 (z= 0.002336). The 
wavelength plotted is in the observed reference frame. 
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Fig. 6. — This figure shows the Fe XXII lines, which are modeled 
with Gaussian profiles. These lines were used as density diag- 
nostics, constraining the density to be rt e < 7 X 10 12 cm - 3 .The 
wavelength plotted is in the observed reference frame, while the 
line identifications are in the rest frame of the galaxy. 
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Fig. 7. — This figure shows how the ratio of Fe XXI lines 
12.38/12.285 change with increasing density. This plot was cre- 
ated using the CHIANTI populations at a temperature of 5 X 10 5 
K. 
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Fig. 8. — The above plot shows the fluorescence Fe Ka line in 
NGC 4051, as well as two narrow Fe K a and /3 lines. The narrow 
lines likely originate in distant material, while the broad fluores- 
cence line likely originates in the inner-most parts of the accretion 
disk. The red line denotes the relativistic fluorescence line char- 
acterized with the Kerrdisk line profile. The blue line does not 
include this line. The narrow Fe K a and /3 features are modeled 
with Gaussian lines. The line markers are drawn in the rest frame 
of NGC 4051. The wavelength plotted is in the observed reference 
frame. 
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Fig. 9. — This is the O VIII Lyman a fluorescence line modeled 
with a Kerrdisk component to describe any relativistic O VIII 
emission. The blue line is model with no relativistic line but in- 
cludes the XSTAR model. The green line includes a relativistic 
line with all of its parameters to vary ( spin, inclination angle, in- 
ner radius, and normalization). The red line denotes the model in 
which the O VIII parameters were tied to that of the Fe Ka line. 
The line is suggested at a 2 a level if the parameters are free to 
vary or tied to that of the Fe Ka line. The wavelength plotted is 
in the observed reference frame, and the data were rebinned with 
a minimum of 25 counts per bin for clarity. 
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Fig. 10. — This figure shows the HETG spectra and the best-fit XSTAR, model components. The parameters for this model are given in 
Table(4] and the model has \ 2 /v = 8177/7516. The absorption features in the residuals as given in Figure[T]are well characterized by the 
XSTAR model, especially in the ranges of 12-18A. The wavelength plotted is in the observed reference frame. 




Fig. 11. — The above figure plots the HETG spectra and the best-fit Cloudy model. The parameters for this model are given in Table[5] 
and the model has x 2 /v = 8713/7516. The wavelength plotted is in the observed reference frame. 



